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ABSTRACT 

The in situ measurement of Comet IP/ Halley and the 'Stardust' returned samples 
of comet Wild 2 showed the presence of a mixture of compact and aggregate par- 
ticles, with composition of both silicates and organic refractory in cometary dust. 
Recently, the result obtained from 'Stardust' mission suggests that the overall ratio of 
compact to aggregate particles is 65:35 (or 13:7) for Comet 81P/Wild 2 (Burchell et 
al. 2008, Meteoritics & Planetary Science, 43, 23). In the present work, we propose a 
model which considers cometary dust as a mixture of compact and aggregate particles, 
with composition of silicate and organic. We consider compact particles as spheroidal 
particles and aggregates as ballistic cluster-cluster aggregate (BCCA) and ballistic ag- 
glomeration with two migrations (BAM2) aggregate with some size distribution. The 
mixing ratio of compact to aggregate particles is taken to be 13:7. For modeling Comet 
IP/ Halley, the power-law size distribution n(a) ~ a~^'^, obtained from re-analysis 
of the Giotto spacecraft data, for both compact and aggregate particles is taken. We 
take a mixture of BAM2 and BCCA aggregates with a lower cutoff size around 0.20/iTO 
and an upper cutoff of about 1/irn. We also take a mixture of prolate, spherical and 
oblate compact particles with axial ratio (E) from 0.8 to 1.2 where a lower cutoff size 
around O.l/zw and an upper cutoff of about 10/im are taken. Using T-matrix code 
for polydisperse spheroids (0.1/xm ^ a ^ lO/im) and Superposition T-matrix code 
for aggregates (0.2/im ^ at, ^ 1/im), the average simulated polarization curves are 
generated which can best fit the observed polarization data at the four wavelengths A 
= 0.365^m, 0.485/im, 0.670//m and 0.684/im. The suitable mixing percentage of ag- 
gregates emerging out from the present modeling corresponds to 50% BAM2 and 50% 
BCCA particles and silicate to organic mixing percentage corresponds to 78% silicate 
and 22% organic in terms of volume. The present model successfully reproduces the 
observed polarization data, especially the negative branch, for comet IP/Halley at the 
above four wavelengths, more effectively as compared to other work done in the past. 
It is found that among the aggregates, the BAM2 aggregate plays a major role, in 
deciding the cross-over angle and depth of negative polarization branch. 
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1 INTRODUCTION 

The study of polarization of the scattered radiation from 
comets, over various scattering angles and wavelengths, gives 
valuable information about the nature of cometary dust. The 
analysis of polarization data gives information about the 
physical properties of the cometary dust, which include size 
distribution, shape and complex refractive indices. 

The in situ dust measurement of Comet IP/Halley 
gave the first direct evidence of grain mass distribution 
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(Mazets et al. 1986). Mukai, Mukai & Kikuchi (1987) and 
Sen et al. (1991) analyzed the polarization data of Comet 
IP/Halley using power law dust distribution (Mazets et al. 
1986) and using Mie theory derived a set of refractive indices 
of cometary grains. The dust distribution function derived 
by Mazets et al. (1986) is actually based on only Vega 2 re- 
sults, while Lamy, Griin & Perrin (1987) derived the grain 
size distribution function for Comet IP/Halley by compar- 
ing the data from spacecrafts Vega 1, Vega 2 and Giotto. 
Much later this dust distribution function was used by Das, 
Sen & Kaul (2004) to analyze the polarization data of a 
number of comets including Comet IP/Halley. 
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Several investigators made useful polarirnetric measure- 
ments of Comet IP/Hallcy through International IP/Halley 
Watch (IHW) filters (Bastien, Menard & Nadeau 1986; 
Kikuchi et al. 1987; Le Borgne, Leroy & Arnaud 1987; Sen 
et al. 1991; Chcrnova, Kisolcv & Jockors 1993). The polar- 
ization data of Comet IP/Halley were analyzed by several 
investigators using Mie theory which assumes the dust par- 
ticles to be spherical (Mukai et al. 1987; Sen et al. 1991; 
Das et al. 2004). However, the naturally occurring cometary 
grains cannot be ideal compact spheres, as required by Mie 
theory. The Mie theory w£is used, as it is more convenient 
and direct, with fewer numbers of free parameters required 
for modeling. Das & Sen (2006) studied the non-spherical 
dust grain characteristics of Comet Levy 1990XX using the 
T-matrix theory. They found that compact prolate grains 
(with axial ratio = 0.486) as compared to spherical grains 
can better explain the observed linear polarization data. As- 
suming an individual cometary grain to be an aggregate of 
several monomers. Das et al. (2008a) again analyzed the ob- 
served polarization data of Comet C/1990 Kl Levy and suc- 
cessfully reproduced the polarization curve through simula- 
tions, where the fit was still better. The X^jn value for the 
aggregates was found to bo 4.2 whereas the value obtained 
by Das & Sen (2006) for compact prolate grains was 5.22. 
Thus it was concluded that aggregate particles can produce 
a still better fit to the observed data as compared to com- 
pact prolate grains. Again, Das et al. (2008b) successfully 
explained the polarization characteristics of comet C/1995 
Ol Hale-Bopp at A = 0.485 ^im and 0.684 /nm using aggre- 
gate dust model. However aggregate dust model used in the 
previous work was restricted to single size of monomer with 
same size parameter at different wavelengths. More recently, 
Das et al. (2010) included the size distribution for aggregates 
and studied the observed polarization data of comet C/1996 
B2 Hyakutake at A = 0.365/im, 0.485/im and 0.684/im. 

It is now well accepted from in situ measurement of 
comets and 'Stardust' returned samples of comet Wild 2 
that cometary dust consists of a mixture of compact par- 
ticles and aggregates (Lamy et al. 1987; Fomenkova et al. 
1999; Horz et al. 2006; Zolensky et al. 2006, Burchell et al. 
2008 etc.). Lasuc et al. (2009) studied comet IP/Halley and 
comet C/1995 Ol Hale-Bopp using a mixture of fluffy aggre- 
gates and compact solid grains. They developed a model of 
light scattering by a size distribution of aggregates of up to 
256 submicron-sized grains (spherical or spheroidal) mixed 
with single spheroidal particles. A good fit of the positive po- 
larization observations of IP/Halley had been obtained by 
them with a power law size distribution (a~^'* with a lower 
cutoff of 0.26 iim and an upper cutoff of 38 jim) with a mix- 
ture of silicates (between 40% and 67% in volume) and more 
absorbing organic material (between 33% and 60% in vol- 
ume) . The fits deduced from their model show that the neg- 
ative polarization branch is not deep enough to match the 
observed polarization data, especially for comet IP/Hallcy. 
Although the fits are found to be good for the positive part 
of the polarization. Recently, Kolokolova & Kimura (2010) 
modeled cometary dust as a mixture of compact particles 
(made of silicate) and aggregates (made mainly of organics 
and IP/Halley like composition). Using a size distribution 
function for compact particles and 256 number of BCCA 
aggregates, they reproduced the polarirnetric data, including 
negative polarization at small phase angles and the positive 



polarization with the maximum value less than 30% at the 
phase angle around 90° and red polarirnetric color. However, 
their model reproduced feature common to 'dusty comets' 
polarization curves but they did not use a chi square fitting 
procedure to compare with observed data for a given comet. 

In the present work, a model for cometary dust with a 
mixture of compact spheroidal particles and aggregates with 
size distribution are proposed to study the observed polar- 
ization data of Comet IP/Halley at A = 0.365/Lim, 0.485/im, 
0.670//m and 0.684;tim. 



2 DUST MODEL 

The in situ measurement of comet IP/Halley and the 'Star- 
dust' returned samples of comet Wild 2 showed the presence 
of a mixture of compact and aggregate particles with com- 
position of silicates and organic refractory. Moreno et al. 
(2007) conducted a systematic approach to test whether a 
collection of compact particles can reproduce the observed 
properties of cometary dust. Using a model of spheroidal 
particles, they found that the axial ratio should be either 
E = 2 (oblate) or E = 0.5 (prolate). The refractive indices 
lie within a range n = 1.6 - 1.7 and k = 0.05 - 0.1. They 
also studied a more complex model based on size distribu- 
tions of irregularly shaped particles composed by a varying 
number of cubes as elementary units. The models consider- 
ing irregularly shaped and compact particles with different 
structures showed results close to the observations. However, 
the weakness of their model of compact structures was that 
the maxima in the linear polarization values did not take 
place in the 90° — 100° phase angle region as observed. Re- 
cently, Kolokolva & Kimura (2010) modeled cometary dust 
as a mixture of compact spheroidal and aggregate particles. 
The compact particles which they considered to be a mix- 
ture of oblate and prolate spheroids with axial ratio within 
the range 1 - 2.5 and aggregates were taken to be Ballistic 
Cluster-Cluster Aggregate (BCCA). 

In the present work, we propose a model which considers 
cometary dust as a mixture of compact and aggregate par- 
ticles. Since the m situ analysis of dust samples exhibits the 
overall ratio of compact to aggregate particles to be 65:35 
(Burchell et al. 2008) , so we take the same value in our analy- 
sis. For modeling comet IP/Halley, we will use a power-law 
size distribution, n(r) = dn/da ~ a~^'^ for both compact 
particles and aggregates, obtained from a re-analysis of the 
Giotto data by Fulle et al. (2000). 

We consider compact spheroidal particles with a size 
distribution from 0.1/im to lO^m. The particles are pre- 
sented by multishaped, polydisperse mixture of spheroids. 
We consider a mixture of prolate, spherical and oblate com- 
pact particles with axial ratio (E) from 0.8 to 1.2. Com- 
putations of light scattering by plain and coated parti- 
cles are made through codes adapted from T-matrix code 
(Mishchenko & Travis, 1996). 

We build the aggregates using ballistic aggregation pro- 
cedure (Meakin 1983, 1984). Two different models of cluster 
growth are taken: first via single-particle aggregation and 
then through cluster-cluster aggregation. These aggregates 
are built by random hitting and sticking particles together. 
The first one is called Ballistic Particle-Cluster Aggregate 
(BPCA) when the procedure allows only single particles to 



join the cluster of particles. If the procedure allows clusters 
of particles to stick together, the aggregate is called Ballis- 
tic Cluster-Cluster Aggregate (BCCA). Actually, the BPCA 
clusters are more compact than BCCA clusters (Mukai et 
al. 1992). The porosity of BPCA and BCCA particles of 128 
monomers has the values 0.90 and 0.94, respectively and 
the fractal dimension of BPCA and BCCA is £) « 3 and 
2, respectively. A systematic explanation on dust aggregate 
model is already discussed in our previous work (Das et al. 
2008a). 

Recently, Shen et al. (2008) considered three different 
classes of clusters distinguished by aggregation rules. These 
are BA (ballistic agglomeration), BAMl (ballistic agglom- 
eration with one migration) and BAM2 (ballistic agglomer- 
ation with two migrations). They developed a set of param- 
eters to characterize the irregular structure of these aggre- 
gates. Actually BA cluster is identical with BPCA cluster. 
The geometry of BAMl and BAM2 clusters are random but 
less porous than BA clusters. The effective porosity (P) in- 
creases from BAM2 BAMl BA. The porosity of BAM2 
structure having 64 number of monomers have the value 
P « 0.5 and the fractal dimension is £) « 3 (Shen et al. 
2008). The aggregates arc taken from web 1 (sec reference). 

In our model, we take same cloud of particles, i.e., same 
type of particles (compact spheroidal and porous (BCCA 
-|- BAM2) particles) and same size distribution (n(r) = 
dn/da ~ o~^'®), to fit the observed data at all wavelengths. 

In our simulation, we divide the present work into two 
phases: 

(i) We first take BCCA aggregates and then mix with 
compact spheroidal particles with 65:35 mixing ratio. The 
result obtained from this modeling will be discussed in Sec- 
tion 4- 

(ii) We then consider more compact aggregate BAM2 
(having porosity (P) ~ 0.50 approximately) which is mixed 
with highly porous BCCA clusters (P ~ 0.9) with some vari- 
able mixing ratio (/3). Then the aggregate mixture is mixed 
with compact particles with 65:35 mixing ratio. Here we 
take composition of both silicate and organic with variable 
mixing ratio (7). 

The free parameters used in the model arc as follows: 

• the mixing ratio (/3) between BCCA and BAM2. 

• the mixing ratio (7) between silicate and organic. 

We use - minimization technique to evaluate the best 
fit values of the above free parameters by the following equa- 
tion: 

^poi-Z^I e;i^) W 

Here, Pobs{6i, A) is the degree of linear polarization ob- 
served at scattering angle 6i (i = 1,2,....,N) and wavelength 
A, Ptheo{di, A) is the polarization values obtained from model 
calculations and Ep{9i, A) is the error in the observed polar- 
ization at scattering angle 9i and wavelength (A). It is also 
observed that this technique of minimization of is quite 
unique. The value of X^in Si^^s the confidence level on our 
best fit values of /3 and 7 and also in the overall fitting pro- 
cedure. Some preliminary work on combined dust model has 
been already reported in Das & Sen (2011). 
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3 COMPOSITION 



The in situ observation of comets, laboratory analysis of 
samples of IDP and remote infrared spectroscopic study 
of comets give useful information about the composition 
of comctary dust. The in situ measurement, of impact- 
ionization mass spectra of Comet IP/Halley's dust, has 
suggested that the dust consists of magnesium-rich sil- 
icates, carbonaceous materials, and iron-bearing sulfides 
(Jessberger et al. 1988; Jessberger 1999). These materials 
are also known to be the major constituents of IDPs (Brown- 
lee et al. 1980). The studies of comets and IDPs have shown 
the presence of amorphous and crystalline silicate minerals 
(e.g. forsterite, enstatite) and organic materials (Hanner & 
Bradley 2004). Laboratory studies have shown that major- 
ity of the collected IDPs fall into one of the three spectral 
classes. These observed profiles indicate the presence of py- 
roxene, olivine and layer lattice silicates. This is in good 
agreement with results obtained from Giotto and Vega mass 
spectrometer observations of Comet IP/Halley (Lamy et al. 
1987). The infrared (IR) measurement of comets has also 
provided important information on the silicate compositions 
in cometary dust. The spectroscopic studies of silicates have 
shown the predominance of both crystalline and amorphous 
silicates consisting of pyroxene or olivine grains (Wooden et 
al. 1999; Hayward, Hanner & Sckanina 2000, Bockelee - Mor- 
van et al. 2002 etc.). Mg-rich crystals are also found within 
IDPs and are predicted by comparing the IR spectral fea- 
tures of Comet C/1995 01 Hale-Bopp with synthetic spectra 
obtained from laboratory studies (Hanner 1999; Wooden et 
al. 1999, 2000). 'Stardust' samples have also confirmed a va- 
riety of olivine and pyroxene silicates in Comet 81P/Wild 2 
(Zolensky et al. 2006). 



It is to be noted that though we used only two free 
parameters j3 and 7 in our model, the refractive indices of 
silicate and organic can be used as other free parameters. 
However, we have limited ourselves to the value taken from 
standard references, because with many free parameters the 
computational time becomes very long. In our computation, 
we take the refractive indices of silicate (especially amor- 
phous pyroxene) from Dorschner et al. (1995). The refractive 
indices of the amorphous pyroxene (Mga^Fei-xSiOs, where 

X is the Mg number, x = x= 0.4, 0.5, 0.6, 0.7, 0.8, 

0.95 and 1.0) are reported by them for different values of x. 
We select x = 0.5 to consider the equal number of Mg and 
Fe in the pyroxene formula. However, we do not claim that 
the choice of a; = 0.5 is unique. The values are calculated 
by linearly interpolating the data obtained from laboratory 
studies. The refractive indices arc given by (1.722,0.101) 
at 0.365/im, (1.692,0.0492) at 0.485^771, (1.673,0.0198) at 
0.670/im and (1.672,0.0185) at 0.684/im. The refractive in- 
dices of organic arc taken from Jenniskens (1993) and the 
values are given by (1.679,0.536) at 0.365/im, (1.842,0.459) 
at 0.485^im, (1.942,0.357) at 0.684^im and (1.949,0.349) at 
0.684/im. The refractive indices of silicate and organic at 
0.485/xm have been already used by Das & Sen (2011) to 
model the optical polarization of comets. 
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4 NUMERICAL SIMULATION 

For modeling comet IP/Halley, we use a power-law size dis- 
tribution, n{a) = dn/da ~ a~^'^ for both compact particles 
and aggregates, obtained from a re-analysis of the Giotto 
data by Fulle et al. (2000). The observed linear polarization 
data of Comet IP/Hallcy is taken from Bastion ct al. (1986), 
Gural'Chuk et al. (1987), Kikuchi et al. (1987), Le Borgne 
et al. (1987), Sen et al. (1991) and Chernova et al. (1993) at 
A = 0.365^im, 0.485^im, 0.670/im and 0.684/im. 

We calculate the scattering properties of spheroidal 
compact particles using T-matrix code (Mishchenko & 
Travis, 1996) for 0.1 < a < 10/jm, where a is the equal 
volume sphere radius of the particle. The step size used to 
integrate the size distribution is 0.01/im. We also calculate 
the scattering properties of the BCCA and BAM2 clusters 
using superposition t-matrix code, which gives rigorous so- 
lutions for ensembles of spheres (Mackowski & Mishchenko 
1996). 

The size of the individual monomer in a cluster plays an 
important role in scattering calculations. These have been 
confirmed by the results of previous work on dust aggregate 
model (Kimura ct al. 2006; Pctrova et al. 2004; Hadam- 
cik et al. 2006; Bertini ct al. 2007; Das et al. 2008a). The 
radius of an aggregate particle can be described by the ra- 
dius of a sphere of equal volume given by av = UmN^^^, 
where N is the number of monomers in the aggregate. In 
the present work, BCCA with 128 monomers and BAM2 
with 64 monomers are taken. As we had computational lim- 
itation and BAM2 cluster takes longer computational time 
compared to BCCA particles, we had to restrict ourselves 
the number of monomers to 64 only for BAM2 particles. In 
our calculation, averages of three random realizations are 
taken for both BCCA and BAM2. The size range of the 
monomer is taken in the range 0.05/im ^ am ^ 0.20/im. 
Thus the lower cutoff radius of the cluster is 0.2 fim and the 
upper cutoff is 1 fj,m. It is to be noted that since the number 
of monomers is fixed in each type of aggregate, the distribu- 
tion in monomer sizes is essentially the size distribution of 
aggregates. For a particular type of aggregate with fixed N, 
the size distribution is just dn/day ~ a~^'^. The step size 
used to integrate the size distribution is O.Ol/im. 

We start calculation considering only BCCA particles 
and then mix with compact spheroidal particles with 13:7 
mixing ratio. It has been checked that the mixing of com- 
pact spheroidal grains and aggregates will not help much in 
producing deeper negative polarization branch beyond 157° . 
It has been also observed from Lasue et al. (2009) that the 
fits deduced from their modeling do not show deep negative 
polarization branch for comet IP/Halley. 

Using aggregate dust models with BAM2 geometry and 
moderate porosity (PwO.6), Shen et al. (2009) reproduced 
albedo and polarization for cometary dust, including nega- 
tive polarization observed at scattering angles beyond 160°. 
To study the effect of BAM2 structure, we now start com- 
putation at A = 0.485/itm with BCCA and BAM2 particles 
with different mixing ratio {f5) and then finally mix with 
compact spheroidal particles having mixing ratio 65:35 (or 
13:7), where 7 is taken to be 3:1. In Fig.l, the polarization 
curves are generated for /3 = 1:3 and 1:1, which actually cor- 
respond to (25% BCCA + 75% BAM2) and (50% BCCA + 
50% BAM2) particles for a size distribution n(a) ~ a~^'^. 



The size range for the aggregates (BCCA and BAM2) and 
compact spheroidal particles is taken to be 0.2 ^ a ^ l.O/Ltm 
and 0.1 < a < 10/xm respectively. We also generate the po- 
larization curves separately with BCCA and BAM2 parti- 
cles. Fig. 1(a) shows the average polarization curve obtained 
from the mixing of compact spheroidal particles and aggre- 
gates (BCCA and BAM2). The mixing ratio between BCCA 
and BAM2 is 1:3. In Fig.l(b), the curve 1 corresponds to 
average polarization curve in the range 150° - 180° obtained 
from the mixing of compact and BAM2 particles only, curve 
2 with /3 = 1 : 3, curve 3 with (3 = 1:1 and curve 4 obtained 
from the mixing of compact spheroidal particles and BCCA 
particles only. 

It is clear from Fig. 1 that the existence of BAM2 struc- 
ture (which is more compact than BCCA) becomes crucial 
in producing the deeper negative polarization branch. Thus 
the introduction of BAM2 aggregate in the aggregate mix- 
ture will help to reproduce the negative polarization well 
which was not possible in previous study by several inves- 
tigators for comet IP/Halley. Actually, interplanetary dust 
particles (IDP) may contain both porous and compact ag- 
gregates. So it will be more realistic if we consider aggregates 
to be a mixture of more compact BAM2 (P ~ 0.5) and more 
porous BCCA (P ~ 0.9) clusters with some mixing ratio /3. 

We now use minimization technique to evaluate the 
best-fitting values of (5 and 7 which can fit to the observed 
polarization data. We have already used this minimization 
technique to fit the observed linear polarization data of some 
comets (Das et al. 2008a,b; Das et al. 2010 and Paul et 
al. 2010), with aggregate models of dust. We need to fine- 
tune the firee parameters and 7 in the model to make the 
best fit to the observed linear polarization data of Comet 
IP/Halley. Some preliminary work on combined dust model 
has been reported by Das & Sen (2011) where they used the 
same technique to simulate the observed polarization data 
of comet Halley at 0A85fim. However their work is limited 
to single wavelength only. 

The best fit values of /3 and 7 are found to be 1:1 and 
78:22 at A = 0.365/xm, 0.485/im, 0.670/im and 0.684/xm. The 
^min '^^1^^^ emerging out from the present analysis are 14.8, 
47.2 and 32.5 at A = 0.365/im, 0.670/im and 0.684/im re- 
spectively, whereas the value obtained by Das & Sen (2011) 
for A = 0.485/im is 56.7. The best-fitting average polariza- 
tion curves at four wavelengths are shown in Fig. 2, Fig. 
3, Fig. 4 and Fig. 5. 



5 DISCUSSION 

The in situ measurement, of impact-ionization mass spec- 
tra of Comet IP/Halley's dust, has suggested that the dust 
consists of magnesium-rich silicates, carbonaceous materials, 
and iron-bearing sulfides (Jessberger et al. 1988; Jessberger 
1999). In our modeling we consider cometary dust as a mix- 
ture of compact and porous particles with composition of 
silicates and organic refractory. The silicate to organic ratio 
coming out from our present work is 39:11 or 78% silicate 
and 22% organic in volume. Thus it can be concluded that 
the silicate composition is dominating in comet IP/Halley 
as compared to organic refractory. 

The negative polarization feature of comet is one of the 
important feature observed in comets. Many comets show 
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Figure 1. (a)The average polarization curve obtained from the 
mixture of compact spheroidal particles and aggregates (BCCA 
and BAM2) for a size distribution n(a) ~ o~^-® at A = 0.485 
/im with /3 = 1 : 3, where the mixing ratio between compact 
and aggregates is taken to be 13:7 and 7 = 3:1. (b) the curve 1 
corresponds to average polarization curve in the range 150" -180° 
obtained from the mixing of compact and BAM2 particles only, 
curve 2 with /3 = 1 : 3, curve 3 with /3 = 1 : 1 and curve 4 
obtained from the mixing of compact spheroidal particles and 
BCCA particles only. 
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Figure 2. Polarization values as observed at wavelength A = 
0.365 /im for Comet IP/Halley by Bastien et al. (1986), Gu- 
ral'Chuk et al.(1987), Kikuchi et al. (1987), Le Borgne et al. 
(1987), Sen et al. (1991) and Chernova et al. (1993). The solid 
curve represents the best-fitting average polarization curve ob- 
tained for compa<;t particles and aggregates (BCCA and BAM2) 
for a size distribution n{a) ~ a~^'^ at A = 0.365 jim. 



negative polarization beyond 157° ((Kikuchi et al. 1987; 
Chernova et al. 1993; Ganesh et al. 1998 etc.). Several in- 
vestigators (Greenbcrg & Hagc 1990; Muinoncn ct al. 1996, 
2007; Tishkovets et al. 2004; Potrova et al. 2004; Hadam- 
cik ct al. 2007 etc.) have discussed the cause of negative 
polarization in comets. Actually, it is important to fit the 
observed polarization data in the positive part as well as 
in the negative branch. Using aggregate dust model. Das et 
al. (2008a,b; Paul et al. 2010) successfully reproduced the 
polarization curves including negative branch observed for 
comets C/1990 Kl Levy, C/1995 Ol Hale-Bopp, C/1996 B2 
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Figure 3. The solid curve represents the best-fitting average po- 
larization curve obtained for compa<;t particles and aggregates 
(BCCA and BAM2) for a size distribution n(a) ~ o~^ ® at 
A = 0.485 iim, taken from Das & Sen (2011). 
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Figure 4. The solid curve represents the best-fitting average po- 
larization curve obtained for compa<;t particles and aggregates 
(BCCA and BAM2) for a size distribution n(a) ~ a~^-^ at 
A = 0.670 fim. 



Hyakutake and C/2001 Q4 NEAT. But it is now well ac- 
cepted that cometary dust consists of compact and porous 
particles. Several investigators studied comet using a mix- 
ture of highly porous aggregates and compact solid grains. 
It has been observed that the plots are showing good fit to 
the positive part of the polarization, but do not show deeper 
negative polarization branch beyond 157". 

In our present work, we take a mixture of aggregates 
(highly and moderately porous) and then mix with compact 
spheroidal grains with 13:7 ratio. It can be noticed from Fig. 
2, Fig. 3, Fig. 4 and Fig. 5 that our modeling can suc- 
cessfully reproduce the positive part as well as the negative 
branch of the polarization at three different wavelengths. 
However, if we just withdraw the BAM2 structure from our 
model, the negative polarization branch will not be repro- 
duced at proper scattering angle values. So it appears that 
the existence of BAM2 particles (which is more compact 
than BCCA) is very important in our grain model as it can 
reproduce a deeper negative polarization branch. Thus our 
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Figure 5. The solid curve represents the best-fitting average po- 
larization curve obtained for compact particles and aggregates 
(BCCA and BAM2) for a size distribution n{a) ~ o~^ ® at 
A = 0.684 ixm. 

modeling can help to explain the polarization characteristics 
of comet IP/Halley successfully at different wavelengths. 

The angular dependence of brightness and linear po- 
larization of compact and porous clusters have been in- 
vestigated by Tishkovots ot al. (2004). They found that 
porous clusters arc brighter almost in the whole angular 
range due to the larger cross section, and these clusters pro- 
duce smoother polarization curves with higher maximum. 
The negative branch in the backscattering direction is shal- 
lower, because the wave interference and near-field effects 
are weaker within the aggregates. However, in compact clus- 
ters, both the interference and the near-field effects play a 
major role in producing negative polarization branch. The 
negative branch is deeper for compact clusters as compared 
to porous clusters. The ininimurn is deeper, and inverted 
angle is shifted to smaller scattering angles. The negative 
polarization is mostly generated by the particles below the 
surface layer of the cluster, where the radiation field is in- 
homogeneous, and the amplitude, phase, and propagation 
direction of the wave change randomly (Tishkovets et al. 
2004). It has been demonstrated by Petrova et al. (2004) 
that the external layer of the clusters plays important role in 
forming the polarization phase curve. The appearance of the 
negative polarization branch and its shape strongly depend 
on the sizes of the scattering elements and on the structure 
of the particle ensemble. In a subsequent work, Shen et al. 
(2009) studied the phase curve and polarization values, as 
produced by BAMl and BAM2 and it was found that more 
compact BAM2 cluster show deeper negative polarization 
branch as compared to BAMl and BCCA. 

Before we conclude we may note that the values re- 
ported earlier for comet C/1990 Kl Levy, showed improve- 
ment of fit in Das et al. (2008a) with aggregate grains (x^ 
value 4.2) as compared to Das & Sen (2006) with compact 
prolate grains (X^ value 5.22). 

It is true that, these two values for comet C/1990 Kl 
Levy, are much lower than the value we obtained in the 
present work for comet Hallcy. It may be noted here that 
in the present work, the fit was made on the data points of 
comet Halley collected from various sources as observed by 
different groups of observers (for example at 0.485 micron, 



we have 86 data points collected from six different groups of 
observers). Such data points collected from diverse groups 
of observers, will always have some inherent scatter in their 
values, as observations are made with different aperture sizes 
and with different sets of filters (with different central wave- 
lengths and FWHM). Besides different groups of observers 
use different instruments, with different spectral responses. 
When we club such data points from various sources, ide- 
ally one should calibrate all the observed data points to take 
into account the above effects. But it is a tedious job and 
normally such corrections are never made as in the present 
case. 

On the other hand, in the two earlier work (Das & Sen 
2006 and Das et al. 2008a) on comet C/1990 Kl Levy, the 
fit was made on data values collected from a single source, 
viz Chernova et al. (1993). Therefore, it is quite natural to 
expect that, the x^ value in the present work on Halley will 
be higher than what has been obtained earlier for C/1990 
Kl Levy. And this is due to the diversities in the sources 
of data points for Hallcy. For example at A= 0.670nm, if 
we exclude the data point (position angle, polarization) = 
(130.0, 12.9), the x^ value just drops from 47.2 to 15.6. 

For comet C/1990 Kl Levy the data points were only 
16 as compared to 86 for Halley at wavelength 0.485(i.m. 
Also for comet Halley we considered a much wider range of 
phase angle values as compared to Levy, which constrained 
our grain model further and increeised the x^ value. What is 
important here to note that, for comet Hallcy no other grain 
model can generate a lower x^ value (indicating a better fit) 
than what has been reported by us in the present work. 



6 CONCLUSION 

(i) A mixture of compact spheroidal grains and aggre- 
gates successfully explains the observed polarization data of 
comet IP/Halley at at A = 0.365/am, 0.485/im, 0.670^m and 
0.684/im. 

(ii) The positive part as well as the negative polarization 
have been successfully generated using the proposed com- 
bined model of cometary dust. 

(iii) With the introduction of distribution of monomer 
sizes and BAM2 cluster (more compact than BCCA), one 
can fit the observed polarization data much better, as com- 
pared to the previous work on comet IP/Hallcy. It is also 
observed that existence of BAM2 structure becomes impor- 
tant in reproducing the deeper negative polarization branch. 

(iv) The best-fitting mixing ratio between BCCA and 
BAM2 {P) is found to be 1:1 (or 50% BAM2 -h 50% 
BCCA). Thus it can be concluded that porous grains in 
comet IP/Halley are composed of both highly and moder- 
ate porous particles. 

(v) The best-fitting mixing ratio between silicate and or- 
ganic particles (7) is found to be 39:11 (or 78% silicate and 
22% organic in volume). 
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